Necessary and sufficient conditions for existence of 
bound states in a central potential 



Fabian Brauf 

Service de Physique Generale et de Physique des Particules Elementaires, 
Groupe de Physique Nucleaire Theorique, Universite de Mons-Hainaut, Mons, 
Belgique 

Abstract. We obtain, using the Birman-Schwinger method, a series of necessary 
conditions for the existence of at least one bound state appUcable to arbitrary 
central potentials in the context of nonrelativistic quantum mechanics. These 
conditions yield a monotonic series of lower limits on the "critical" value of the 
strength of the potential (for which a first bound state appears) which converges to 
the exact critical strength. We also obtain a sufficient condition for the existence 
of bound states in a central monotonic potential which yield an upper limit on 
the critical strength of the potential. 
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1. Introduction 

The problem of finding upper and lower limits on the number of bound states of a 
given potential is become a classical problem since the pioneer works of Jost and Pais 
in 1951 pP and Bargmann in 1952 '7. They obtained, for the first time, a necessary 
condition for the existence of bound states in a central potential which can be obtained 
from the following upper limit on the number of ^-wave bound states (setting Ni to 

1) 

Ni<,^^j^drr\V-{r)\. (1) 

In this inequality, (r) is the negative part of the potential obtained by setting its 
positive part to zero and £ is the angular momentum. Note that we use the standard 
quantum-mechanical units such as ?i = 2m = 1, where m is the mass of the particle. 
This upper limit called Bargmann-Schwinger upper limit in the literature, was 
the starting point of intensive studies and a fairly large number of upper and lower 
limits on the number of bound states for various class of potentials was found, see for 
example ,3. A A <g| i7i ig» .9> M .H. .12. .Igt |14| .l^, ,lg. ,17. .1§. .19. .2Q, ,21, . 

An important theorem for classifying these results was found by Chadan |H] and 
gives the asymptotic behavior of the number of bound states as the strength, g, of the 
central potential goes to infinity: 
(,1/2 /-oo 

Nk^ — / drvirf/'^ as g ^ oo, (2) 

Jo 

where the symbol w means asymptotic equality and V~{r) — —gv{r). This result 
implies that any upper and lower limit from which could yield cogent results should 
behave asymptotically as g^^^. More importantly, the relation Q gives the functional 
of the potential, that is to say the coefhcient in front of g^^^, that appears in the 
asymptotic behavior. The upper limit is proportional to g instead of g^^^ and is 
not very stringent for strong potentials. Upper and lower limits featuring the correct 
3^/^ dependency was first obtained in the Ref. Upper and lower limits featuring 
the correct asymptotic behavior ^ was first derived in Refs. |19[l2l)j . In practice, the 
asymptotic regime is reached very quickly when the strength of the potential is large 
enough to bind two or three bound states. 

The situation is completely different when one consider the transition between 
and 1 bound state and in particular upper and lower limit on the "critical" value of 
the strength of the potential, gc, for which a first bound state appears. In this case, 
there is no theorem to know in advance which limit yield the most stringent restriction 
on gc- It is then of interest to obtain various limits, since the limit yielding the most 
stringent restriction change from one potential to another. 

In Section |21 we obtain a series of necessary conditions for the existence of at 
least one bound state, applicable to arbitrary central potentials, which converges to 
the exact critical strength. In Section |31 we present a sufficient condition for the 
existence of bound states in a central monotonic potential. In Section 0] we perform 
several tests of the cogency of the limits presented in this article and we compare them 
to some previously known results and to the exact results. 
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2. Necessary conditions 

The necessary conditions for the existence of bound states derived in this section is 
obtained with the help of a simple extension of the Birman-Schwinger method. Birman 
P] and Schwinger |^ have shown how to obtain an upper limit on the number of bound 
states once the Green function of the kinetic energy operator of a wave equation is 
known. We recall briefly the main line of the method applied to the radial Schrodinger 
equation for completeness; for more details see the original articles 
The Schrodinger equation for a central potential V{r) reads 

' 'Ke{r) = {E-V{r))u,{r). (3) 



dr 
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The zero-energy Schrodinger equation can be written under the form of the following 
integral equation 

u,{r) = - / dr' gt{r, r') V{r') ue{r'), (4) 







where gi{r,r') is the Green function of the kinetic energy operator and is explicitely 
given by 

9^iry) = ^/+^r-\ (5) 

where r< = min[r, r'] and r> — max[r, r'] . Since the purpose of the method is to obtain 
an upper limit on the number of bound states, we can replace V{r) by — [^'"(r)! where 
V~{r) is the negative part of the potential obtained by setting the positive part of the 
potential equal to zero. Indeed, a decrease of the potential in some region must lower 
the energies of the bound states and therefore cannot lessen their number. Moreover, 
we introduce the parameter < A < 1 by the substitution \\V~{r)\. As 

A increases from 0, we reach a critical value, Ai, at which a bound state first appears 
with a vanishing binding energy, E = 0. With further growth of A, the energy of this 
state decreases until we reach a second critical value, A2, at which a second bound 
state appears and so on. When A has attained the value unity and, Ajv^ < 1 < Xn^+i, 
there are Ne bound states. 

We now introduce, to obtain a symmetrical kernel, a new wave function as 

Mr) = \V-ir)\'/'u,ir). (6) 

The equation (0J becomes 



00 



X~'Mr) = / dr' Ki{r,r')Mr'), (7) 





where Ki{r,r') is given by 

K,{r,r') = \V~ir)\'/^geir,r') \V-ir')\'/^. (8) 

The kernel being positive, we have 0<Ai<A2<---<AAr<l and < Afc < 00 (Afc 
denotes each eigenvalue of 0). It is well known that the trace of the iterated kernels 
equals the sum of the eigenvalues of the integral equation as follow 

drX(")(r,r), (9) 
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with 



where the iterated kernel K^"'\s,t) is given by 

K^;^\s,t) = / du Ki{s,u)K^;''^\u,t), (10) 

Kf\s,t)=K,{s,t), (11) 
and n = 1, 2, . . .. Now it is plain that the following inequalities hold 

where Ni is the number of wave bound states. From and 112() we find 

that an upper limit on the number of £-wave bound states of the Schrodinger equation 
is given by 

/■oo 

Ne< drK''l'\r,r). (13) 



In his article, Schwinger consider only the case n = 1 for the equation H13|) which 
yields the Bargmann-Schwinger upper limit Indeed, greater values of n would 
yield upper limits which possess a worse dependency on the strength of the potential 
g than the upper limit and which would be very poor for strong potentials. But 
it appears that, as illustrated in Sectional the larger is the value n the better is the 
lower limit on the critical value of strength of the potential. 

The necessary conditions for the existence of £-wave bound states obtained from 
read respectively for n = 1, 2, 3: 
1 /"^ 

drr\V-{r)\>l, (14) 



2^+1 
2 



(2^+1)2 
6 



poo pri 

/ dr,rY''\V-{ri)\ drs +2 |y-(r2)| > 1, (15) 

/ dr,r-^'\V-{r,)\ dr^r^lV" {r^)] 
Jo Jo 



(21 +1) 

drsrl'+'\V-ir,)\>l, (16) 



The improvements of the lower limits on gc implied by the relations H15|l and 1)161) over 
the lower limit inferred from the well known relation H14|) are illustrated in Section 0] 
for a square-well potential and an exponential potential. 

Let us end this section by noting that the procedure employed here yield also 
a necessary condition for the existence of bound states analogous to the condition 
obtained by Glaser et al. [TU] 

(p - i)p-^r(2p) 
(2^ + i)2p-ippr2(p) 

where p > 1 must be chosen to optimize the result. Indeed, for ^ > 0, we can use the 
n times the Holder inequality in the relation (|13(l and taking n going to infinity (see 
122] for more details) we obtain 

i2i+l)pip-l) TP-i 



drr^P"! |y-(r)|f > 1, (17) 



, , drr^P-^\V-(r)\P > I. (18) 

[p^£ + 2)ie-l) + 3p-l\ Jo I Wl - K ) 

The constant in front of the integral is unfortunately always greater than the constant 

appearing in the necessary condition H17|) . and the relation (I18|l is thus always less 

stringent. 
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3. Sufficient condition 

The sufficient condition is obtained with the help of a generaHzation of the comparison 
theorem proved recently and where the comparison potentials intersect (Theorem 7 of 
Ref. 23 ). The new theorem reads 

Theorem. // two monotonic potentials Vi(r) and V2{r) cross twice for r > at 
T — ri,r2 (ri < r2) with 

(i) Vi{r) < V2{r) for < r < ri and 

(li) r dy[Vi{v)-V2{v)W <Q, 
Jo 

then El < E2, where Ei^2 ctfe the ground states of the potentials Vifl{r). 
As comparison potential V2{r), we choose a simple square- well 

V2{r) = -V^e{R-r), (19) 

where 9{x) is the Heaviside function. Moreover, we choose this potential such as a 
zero-energy bound state exists: VqR^ — 7r^/4. This implies that the potential Vi{r) 
possesses at least one bound state. For this particular choice of V2{r) we have r2 = R- 
We write the potential Vi (r) under the form 

Vi{r)^^gs-^v{rls,k), (20) 

where k are the other parameters of the potential. The hypothesis {ii) above yields 
the following upper bound 5"^ on the critical coupling constant gc 

u 

12 dyy^v{y,k)' 

where a = R/s. The best restriction is obviously obtained with the value of a 
minimizing the right-hand side of l|21|) . The upper limit can thus be written as 

3c ^ = 2 ^ .V (22) 
12 v(a, k) 

where a is the unique solution of 

dyy^v{y,k)=aMa,k). (23) 



The definition H23(l of a has a simple geometric signification which allows to remark 
that a > max[y^?;(j/, k)]. 

Obviously, we have used a very particular comparison potential V2{r) to write 
a neat formula for the upper limit on the critical coupling constant gc- In practice, 
a better upper limit could be obtained by the use of a more appropriate comparison 
potential for which the exact value of the critical coupling constant is known (and for 
which the conditions (i) and (ii) apply!). 

4. Tests 

The first potential we consider to test the limits presented in the previous sections is 
a square-well potential that we write in the convenient form 

Vir) = -gR-^0{l-r/R). (24) 
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Table 1. Comparison between the exact values of the critical coupling constant 
Qc of a square well potential for various value of £ and the lower limits, g]P < gc, 
obtained with the relations 1251 - 1^71 . the lower limit obtained with the relation 
llijj with n = 4 and Ni = 1 (calculated numerically) and the lower limit obtained 
with the formula 1171 (with the optimal value of p). 



e 


n = 1 


n = 2 


n — 3 


n = 4 


Eq. fT7|l 


Exact 





2 


2.4495 


2.4662 


2.4672 


2.3593 


2.4674 


1 


6 


9.4868 


9.8132 


9.8592 


9.1220 


9.8696 


2 


10 


18.708 


19.895 


20.120 


18.454 


20.191 


3 


14 


29.699 


32.383 


32.981 


30.245 


33.217 


4 


18 


42.214 


47.064 


48.272 


44.425 


48.831 


5 


22 


56.089 


63.788 


65.868 


60.947 


66.954 



Table 2. Comparison between the exact values of the critical coupling constant 
Qc of an exponential potential for various value of £ and the lower limits, g].° < gc, 
obtained with the relations 1251 - 1^71 . the lower limit obtained with the relation 
I13i with n = 4 and Ni = 1 (calculated numerically) and the lower limit obtained 
with the formula 1171 (with the optimal value of p). 



i 


n = l 


n = 2 


n = 3 


n = 4 


Eq. 117^ 


Exact 





1 


1.4142 


1.4422 


1.4453 


1.4383 


1.4458 


1 


3 


6.2700 


6.8546 


6.9913 


7.0232 


7.0491 


2 


5 


13.145 


15.257 


15.804 


16.277 


16.313 


3 


7 


21.593 


26.265 


27.364 


29.218 


29.259 


4 


9 


31.363 


39.616 


41.296 


45.849 


45.893 


5 


11 


42.297 


55.120 


57.480 


66.173 


66.219 



The sufficient condition H22 |l -()23 |1 . applicable only for £ = 0, is saturated for this 
potential (with a = 1) and thus leads to the exact result. The necessary conditions 
p4|l - (|16|) give the following lower limits 

= 2(2^+1), (25) 
gi° = (2^+l)[2(2^ + 3)]i/2, (26) 
gl° = (2^+ l)[(2^ + 3)(2^ + 5)]^/^ (27) 

The comparison between the new lower limits on the limit (|17|l and the exact 
results is reported in Table ^ sj^id shows that the new limits are quite cogent and 
converge quickly to the exact result especially for small value of £. 

The last test is performed with an exponential potential written as 

V{r) = -gR-^ exp(-r/i?). (28) 

For £ = 0, the sufficient condition (|22|I - H23|) leads to g"? — 2.118 while the exact result 
is given by gc = Zq/A = 1.4458 {zq = 2.4048 is the first zero of the Bessel function 
Jo{x)). The upper limit is not very stringent for this potential because the comparison 
potential that we choose (a square well) is very different from an exponential potential. 
The upper limit yields more cogent result, for example, for a Wood-Saxon potential. 
For an exponential potential a better upper limit can be obtained with the Calogero 
lower bound 5 : g^^ = 1.677. 
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The comparison between the news lower hmits on g^, the Uniit (|17|) and the exact 
result is reported in the Table |21 The new lower limits on are quite cogent and 
converge quickly to the exact results especially for small value of £, but this convergence 
is slower than in the case of a square well potential. 
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